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Abstract. The time evolution of giant {D > \ Mpc) lobe-dominated 
galaxies is analysed on the basis of dynamical evolution of the entire 
FRII-type population. 



1. INTRODUCTION 

One of the general questions concerning the largest radio sources 
is: do they reach their extremal giant sizes due to (i) exceptional 
physical conditions in the intergalactic medium, (ii) extraordinary 
intrinsic properties of the AGN, or simply (iii) because they are ex- 
tremely old? To answer this question, a number of attempts were 
undertaken to recognize properties other than size which differenti- 
ate giants from normal-size sources. As a result of those attempts, a 
recent notion is that the largest sizes result from a combination of 
the above properties. 

In this contribution we analyse whether properties of giant ra- 
dio galaxies observed in a selected representative sample can be ex- 
plained by a model of the dynamical evolution of classical double 
radio sources in cosmic time, and what factor (if there is a one) is 
primarily responsible for the giant size. Two recent analytical models, 
published by Kaiser et al. (1997) [hereafter KDA] and Blundell et al. 
(1999), are very convenient for this purpose. These two models differ 
in predictions of the time evolution of the source luminosity (there is 
no space in this contribution to go into details). In summary, basic 
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physical parameters, i.e. the jet power Qjet, the central density of 
the galaxy nucleus po, the energy density and pressure in the lobes 
or cocoon (uc and pc), and the total energy of the source Etot are 
derived from the model for each member of the sample to fit its age, 
redshift, radio luminosity, projected size, and axial ratio. Next, these 
parameters are compared with (1) the relevant parameters derived 
for normal-size sources in a comparison sample, and (2) the param- 
eters determined from observational data, i.e. the age, equipartition 
energy density Ueq, equipartition energy Ueq, etc., calculated under 
'minimum energy' conditions. 

A description of the observational data used, application of the 
analytical model, fitting procedure, etc. (cf. Machalski et al. 2003), 
are beyond the scope of this contribution. 

2. RESULTS OF THE MODELLING 

Jet power Qjet and core density po 

A distribution of the above 
parameters on the log((5jet)-log(/9o) 
plane is shown in Fig. 1. As both 
parameters should be independent, 
we test whether the observed dis- 
tribution is or is not biased by pos- 
sible selection effects. The data ob- 
tained imply that Qjet correlates, 
in order of the significance level, 
with the luminosity Pi. 4, redshift z 
and age t. Calculating the Pear- 
son partial correlation coefficients, 
we found no significant correlation 
between Qjet and po when z (or P1.4) 
and t are kept constant. 

In view of the above, one can 
see from Fig. 1 that (i) among the 
sources with a comparable Qjet, giant sources have an average po 
smaller than its corresponding value in normal-size sources, (ii) gi- 
ants have at least ten times more powerful jets than much smaller 
low-luminosity sources of a comparable po, and (iii) for a number of 
sources in the sample, the derived values of their fundamental pa- 
rameters Qjet and po are very close, while their ages and axial ratios 



S 38.5 







« CygA ; 























\ 









* 


* 


Giant * ; 


* * 




low 7 i".' 


* 




high / ' 
low P * - 









-23.5 -22.5 

log(p„[kgm"^]) 



Fig. 1. Plot of the jet power Qjet 
against the core po- The 'clans' of 
sources with similar Qjet and po are 
marked by the large circles 
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are evidently different. Thus in view of the model assumptions, they 
may be considered as 'the same' source observed at different epochs 
of its lifetime. Such bunches of three to five sources (hereafter called 
'clans') are indicated in Fig. 1 with the large circles. 
Relation between D, Qjet, Po, t, and (1 + z) 

In view of the dynamical model applied and as a result of the 
Pearson partial- correlation coefficients calculated between those pa- 
rameters, we find that the linear size of a source strongly depends on 
both its age and the jet power, while the correlation with age is the 
strongest. However, the size also anti-correlates with central density 
of the core. That anticorrelation seems to be weaker than the cor- 
relations with Qjet and t and become well pronounced only when all 
three remaining parameters {Qjet, t and z) are kept constant. 

3. EVOLUTIONARY TRACKS OF SOURCES 

In the papers of KDA and Blundell et al. the tracks of radio 
luminosity P„ versus linear size D were derived for imaginary sources 
with assumed values of Qjet, po, z and other free parameters of the 
model. In our approach we are able to calculate such evolutionary 
tracks for actual sources. In Sect. 2, the 'clans' of sources have been 
pointed out. Three of six clans are marked in Fig. 1. Since the dy- 
namical model assumes constant Qjet during a source lifetime, and 
a nucleus density po is a priori constant, members of such a clan can 
be considered as 'the same' source observed at a number of different 
epochs throughout its life. The observed parameters of these mem- 
bers can verify predictions of the model. However, fits of the P-D 
tracks predicted with the original KDA model to the observed pa- 
rameters of sources have appeared unsatisfactory. Much better fits 
of the modelled tracks to observational data of the 'clan' members 
are found with the cocoon's axial ratio (AR) evolving in time. This, 
in turn, implies a time evolution of a ratio of the jet head pressure 
to the cocoon pressure, Phc- 
The tracks logPi.4-logZ) and logUeq-log Etot 

The evolutionary P - D tracks for the three clans are shown in 
Fig. 2a with solid curves. The markers of the same time are put on 
these tracks. The dashed curves show the relevant tracks calculated 
from the original KDA model, i.e. with a constant AR taken as 
the mean of axial ratios in a given clan. It is clearly seen that the 
evolving AR much better fits the observed changes of P on D. 
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Fig. 2(a). Evolutionary P — D tracks and (b) u^q — Etot tracks fitted for three 
clans of sources with evolving axial ratio AR (solid curves). The markers of the 
same predicted age on each curve are connected with dotted lines. The members 
of each clan are marked with different symbols. Their actual age is indicated 
by a number behind the symbol. The dashed curves indicate relevant tracks but 
calculated with a constant AR, as in original KDA model 

The model also allows a prediction of the source's evolution 
on the energy density-total energy plane which are shown in Fig. 
2b. It is worth emphasizing that the predicted evolutionary Uc-Etot 
tracks are steeper and curved in respect to those expected from the 
original KDA model. Moreover, the steepening increases throughout 
the source lifetime. This is caused by the non-constant adiabatic 
losses and inflation of the cocoon in time, as well as a faster decrease 
of the cocoon pressure in very large sources. Quantitatively this 
process is evaluated by a substitution of the evolving (increasing) 
value of the pressure ratio Phc{t) into equation describing Etot (cf. 
Machalski et al. 2003). 
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